Evaporative Cooling in Semiconductor Devices 
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We discuss the theory of cooling electrons in solid-state devices via "evaporative emission." Our 
model is based on filtering electron subbands in a quantum-wire device. When incident electrons 
in a higher-energy subband scatter out of the initial electron distribution, the system equilibrates 
to a different chemical potential and temperature than those of the incident electron distribution. 
We show that this re-equilibration can cause considerable cooling of the system. We discuss how 
the device geometry affects the final electron temperatures, and consider factors relevant to possible 
experiments. We demonstrate that one can therefore substantial electron cooling due to quantum 
effects in a room-temperature device. The resulting cooled electron population could be used for 
photo-detection of optical frequencies corresponding to thermal energies near room temperature. 
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PACS numbers: 73.50.Lw,73.23.-b 



I. INTRODUCTION 



As electronic devices become smaller, they leave the 
regime of classical physics and enter the realm of quan- 
tum physics. Many classical quantities such as resistance 
must be reinterpreted for systems on a mesoscopic scale. 
One such classical concept is that of the refrigerator: a 
device that uses an external source of work to cool a gas. 
In this paper we consider whether this classical concept 
can be applied to an electron gas so that one could cool 
such a gas by applying a voltage to a device. 

There are many ways to cool electrons in a condensed- 
matter system. For example, thermoelectric coolers 
based on the Peltier effect are available commercially. 
A different kind of electron-cooling mechanism in semi- 
conductor devices is based on a quasistatic expansion of 
a two-dimensional electron gasiS Still other possibilities 
include taking advantage of many-body effects that can 
lead to liquid/gas phase transitions in the electron pop- 
ulation in a semiconductor quantum wellA^ 

In this paper we investigate electron cooling in a meso- 
scopic solid-state device using evaporative emission. This 
method entails removal of ( "filtering" ) electrons from a 
high-energy subband of a many-electron system, followed 
by relaxation of the remaining electrons to a temperature 
lower than that of the initial system. Evaporative cool- 
ing is widely used in bosonic systems^ But this method 
is harder to implement for fermionic systems, as we shall 
discuss below. 

In section [H] we describe the theory of a two- 
dimensional device to cool electrons in quantum wires. 
We use the Landauer formula^i^ to analyze the cooling 
properties of these devices. This formula was originally 
developed to explain the transport properties of electrons 
in a quantum device. It relates these properties to the 
quantum mechanical scattering amplitudes for electrons 
that pass through the device. To calculate these ampli- 
tudes, we use an extension of R-matrix theory2ii2*Ai that 
we summarize in the Appendix. In section IIIII we use 
this theory to calculate cooling properties of several two- 
dimensional devices. We begin with a simple T-junction 



device and show that by optimizing its design we can 
achieve electron cooling. We improve upon this result 
by switching to a "plus-junction" design, which can give 
up to 15% cooling. In section HVI we discuss applications 
and realistic parameters for a device to cool electrons, 
and in|7|we summarize our key results and describe fu- 
ture research. 



II. ELECTRON COOLING IN 
TWO-DIMENSIONAL QUANTUM DEVICES 

Our theoretical approach is analogous to the work- 
ing principle of the classical Hilsch vortex tube, 13 which 
uses a T-shaped assembly of pipes to separate high- 
pressure air into a high-temperature system and a low- 
temperature system. This separation does not violate the 
Second Law of Thermodynamics, because the system is 
driven by an external force. 

We use a similar idea to cool electrons in a quan- 
tum mechanical system. The simplest such device uses a 
T-shaped assembly of quantum wires to remove higher- 
energy electrons from an electron gas at fixed tempera- 
ture. Figure^is a schematic of such a configuration that 
defines the regions of the device. We assume that our 
device is formed from a quantum well whose thickness is 
sufficiently small that the device can be considered two- 
dimensional. That is, we assume that the device confines 
electrons to a layer of thickness zq such that the confine- 
ment energy associated with motion in the z-direction is 
much larger than any other energy in the problem. (This 
confinement energy is ft 2 7r 2 /2m* z 2 , where m* is the elec- 
tron effective mass.) 

There are three leads in the T-junction: the input lead, 
output lead, and sidearm, as shown in Fig^ We shall 
label physical quantities by subscripts "i," "o," and "s" 
respectively; for example, we denote the widths of the 
leads by tWj, w , and w s . Electrons are injected into this 
device through the input lead, and in the input region are 
in thermal equilibrium at an initial temperature Ti and 
chemical potential /i,-. Filtering of higher-energy elec- 
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FIG. 1: Regions in a two-dimensional T-junction device. The 
"ballistic region" is the part of the device within a distance 
comparable to the electron-electron scattering length (t e - e ) 
of the junction. The "adiabatic region" is within the electron- 
phonon scattering length (£ e - p h) of the junction. It is in the 
adiabatic region that we achieve cooling. At larger distances, 
in the "equilibrium region," the electrons have returned to 
the temperature of the lattice. 



A. Input electron densities and populations 

Provided the electrons in each lead are in thermody- 
namic equilibrium, we can treat them as an ideal Fermi 
gas. For lead £, therefore, the subband electron density 
(per unit volume) in an open (energetically accessible) 
subband n is 



p M{E,T t ,M) = f(E-,T t ,M)vV>(E) 



(1) 



(£) 

for n such that the subband energy obeys e n < E. In 
this equation, the density of states is 



V^{E)=[E- 



e (Z)1-V2 



(2) 



The quantity E is the dimensionless electron energy mea- 
sured in units of the lowest subband energy of the input 
channel, £\ = ti 2 ir 2 /2m* wf . We choose the zero of en- 
ergy at the energy of the ground transverse state in the 
input lead. We measure all other subband thresholds 
relative to this energy and in units of £ \ , so that 



trons from the initial electron gas occurs in the scattering 
region. The rate of scattering into the sidearm depends 
upon the electron energy: if the subband energy in lead of 
width w is given by E n = h 2 -K 2 n 2 /2mw 2 , then the "force" 
electrons exert on the sides of the lead, F = —dE/dw, 
is larger for larger n. Thus when an electron encounters 
the sidearm, the higher-subband states are more likely to 
squirt down the sidearm. Alternatively, one can see that 
the higher-subband electrons scatter preferentially to the 
sidearm because they are in states with wave functions 
that are linear combinations of plane waves with larger 
transverse momenta. 

Electrons that scatter forward into the output lead pro- 
ceed into the ballistic region, where the electron popula- 
tion is determined entirely by the product of the initial 
electron distribution and the scattering probability. Next 
the electrons enter the adiabatic region, where they ex- 
change energy among themselves and so relax to a tem- 
perature T Q and a chemical potential /i D . In principle, the 
values of these properties can be calculated from conser- 
vation laws for energy and particle number. In section 
IIIII we report such calculations and show that for some 
device geometries the temperature T a is less than the ini- 
tial temperature: this temperature decrease is the desired 
cooling effect. 

Finally, at large distances (as determined by the 
electron-phonon scattering rate) the electrons will return 
to equilibrium with the lattice at the initial tempera- 
ture Ti, this re-equilibration occurs in the equilibrium 
region. To see cooling, one must measure the electron 
temperature before they get to this region 

Below we define our notation and describe how we cal- 
culate the electron distributions in the input and output 
leads. We also define and describe the calculation of a 
quantitative measure of electron cooling in the device. 
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nr — y - 1. 

wi 



(3) 



The occupation probability for lead £ is given by the 
Fermi-Dirac distribution function for total electron en- 
ergy E = E n e) + e£\ 



f{E;T,m) 



1 



e {E~p, e )/k B T _|_ i ' 



(4) 



where fi£ is the electrochemical potential in lead £, En 
is the longitudinal kinetic energy of the electron, and 
ks is Boltzmann's constant. The corresponding subband 
population in lead £ is the integral of the subband den- 
sity (Eq. QJ) over all allowed total electron energies E: 



p${E;T tilH )dE. 



(5) 



The total population in lead £ is the sum of the subband 
populations for that lead, 



JVW(T /)W )=gjvW(T <)M/ ). 



(6) 



B. Transmitted-electron densities and populations 

Immediately upon leaving the scattering region the 
transmitted electrons are in a highly noncquilibrium dis- 
tribution and cannot be characterized by a temperature 
or a chemical potential. By the time these electrons have 
traveled a distance along the lead comparable to several 
times their relaxation length, they have come to equi- 
librium at T a and p , and it is meaningful to describe 
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them by a Fermi-Dirac distribution function f(E; T Q , p a ). 
These output-lead properties refer to the electron popu- 
lation in the adiabatic region, where the electrons are in 
thermodynamic equilibrium. Given and pi for the in- 
put lead, our goal is to determine the values of T D and p a 
for the output lead that yields the lowest T < T; i.e., 
which maximizes cooling of transmitted electrons. 

For a given T, and pi, we can determine the tempera- 
ture in the output lead T a by requiring that the number 
of electrons in the output lead at equilibrium equals the 
number of electrons transmitted into this lead (conserva- 
tion of electrons in the output lead) . To set up equations 
to implement this strategy, we must define subband and 
lead populations in terms of electrons transmitted from 
the input lead into the output lead. The (state-to-state) 
density of electrons transmitted from an open subband n 
of the input lead i into an open subband n! of a lead £' is 

p%{E;T ulH ) = $(E,T u h)T** (7) 

for E > e™, a *, where T i' 1 is the transmission coefficient 
from subband n in lead i to subband n' in lead £'. The 
restriction that the total electron energy E be greater 
than or equal to 

e^=max{e»e^} (8) 

ensures that subbands n and n' are both open; were 
this restriction violated, then the transmission coefficient 
T n , ' l n would be undefined. 

The population of electrons transmitted into sub- 
band n' of the output lead, the transmitted subband den- 
sity, is 

oo />oo 

N°j{T i , l i i ) = Y J / P y >n (E;T uN )dE. (9) 

Hence the total population of electrons transmitted into 
the output lead, the transmitted density, is 

oo 

N°' i (T i , tH )=Y / N^(T i , f i i ). (10) 

n'=l 



C. The cooling parameter 

For a given Ti and pi, we can determine the equilib- 
rium temperature in the output lead T Q by solving simul- 
taneously the equation for conservation of the number of 
electrons, 

N^ \T ,p )=N ' i {T i ,p i ), (11) 



and the equation for conservation of electron energy in 
this lead: 

(E)M(T ,p ) = (E}°> i (T i ,p i ). (12) 

In Eq. (|12|) the energy of the electrons at equilibrium in 
the output lead is 

°° POO 

(E)^(T ,fx )=J2 Ep^)(E;T ,p )dE, (13) 

n'=l *C' 

where the subband density in the output lead, 
p { °) {E;T Q , p Q ) is defined by Eq. (TJ. The energy of the 
electrons transmitted into the output lead is 

oo oo roc 

(E)°> i (T i ,p i )=J2Yl / Ep°r;UE;,Ti,iH)dE, 

n'=lt»=l ^ e '"' a ^ 

(14) 

where the transmitted-electron density p£ ' n (E; , Ti, pi) is 
given by Eq. Q. Using Eqs. (JTDJl, and (JT4J we calculate 
the number of transmitted electrons and their energy. We 
then calculate the equilibrium temperature and chemical 
potential that would give the same total number and en- 
ergy. This calculation produces the parameters T a and 
p for the output lead. 

As a measure of the effectiveness of a given device for 
cooling electrons, we define the cooling parameter 

fm \ — T (Ti, pi) , 1C , 
V{Ti,pi) = = . (15) 

^ i 

If r\ > 1, the device heats electrons. Our goal, therefore, 
is to determine the device geometry and initial electron 
properties Ti and pi that minimize rj < 1. 



III. RESULTS 

To achieve optimum cooling the higher-subband elec- 
trons should scatter into the sidearm, and the lower- 
subband electrons should scatter into the output lead. In 
a real device not all higher-subband electrons will scatter 
into the sidearm, and not all the lower-subband electrons 
will scatter into the output lead. The probabilities for 
electron scattering, as quantified in transmission coeffi- 
cients, depend on system properties such as the geome- 
try, and on scattering potentials. We consider "perfect 
devices" that have no impurities, so electrons are scat- 
tered only by the boundaries of the device. We further 
assume that the potential energy of the electrons in the 
leads is zero; this simplifying assumption is not essential 
to either the cooling effect or to our formalism. We can 
alter the scattering of electrons by changing the ratio of 
the width of the sidearm to that of the input lead (Fig.0 
or by changing the geometry altogether. 
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A. T-junction cooling devices 

To calculate the cooling parameter rj we need to know 
the population of electrons in each subband. To deter- 
mine this quantity we must know transmission coeffi- 
cients from a state in the input lead to states in the out- 
put lead. To calculate these transmission coefficients we 
use a generalization of R-matrix theory that we summa- 
rize in the Appendix. To determine the results reported 
here, we calculated cooling parameters 77 using transmis- 
sion coefficients in a T-junction device for various ratios 
of the width w s of the sidearm to the width Wi of the in- 
put lead, keeping the width w a of the output lead equal 

to Wi. 



1. A T-junction device with w s = Wi = 1.0 

We first consider a T-junction device in which all leads 
have the same width: Wi = w s — w Q . Figure Ufa) 
shows "state-to-lead" transmission coefficients for scat- 
tering into the output lead of electrons in different sub- 
bands of the input lead. These coefficients are sums 
over all energetically accessible ("open") subbands n 
of the output lead of state-to-state transmission coeffi- 
cients 7^f ' l n . [see Eq. [7] from a given state of the in- 
put lead. This figure illustrates the loss of some higher- 
subband electrons from the initial electron distribution. 

The cooling parameter 77 for this case is shown in 
Fig. [2Jb) for different values of initial temperature Ti 
with the initial chemical potential fii — 0. (Note that we 
measure the energy in terms of £{, and all the energies 
are measured from £ \ . So fii — means that the external 
potential of the system is such that the Fermi energy is 
Ep = H 2 TT 2 /2m* wf .) For this geometry r\ > 1 for all 
initial temperatures the cooling parameter. That is, this 
device heats electrons — the opposite of the desired effect. 

This case is important because it demonstrates that 
even if high-energy electrons are lost due to scattering, a 
compensatory loss of low-energy electrons may produce 
an overall heating effect. Loss of low-energy electrons 
opens gaps in the electron distribution at low energies. 
Higher-energy electrons can then relax into these newly 
accessible low-energy states, with the resulting energy 
difference liberated as thermal energy. If this happens, 
then the loss of low-energy electrons will heat the system. 
Even a small dip in the scattered-electron distribution at 
low energies will significantly affect the final tempera- 
ture. To cool electrons, therefore, it is not sufficient to 
merely scatter higher-energy electrons. We must scatter 
thermally excited electrons but not significantly scatter 
electrons in lower-energy subbands. 



2. Alternative T-junction geometries 

To determine whether a T-junction device can cool 
electrons, we now consider several widths u, of the 
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FIG. 2: Upper panel: State-to-output-lead transmission co- 
efficients for electrons in a T-junction device with un = w s = 
w . The curves correspond to different subbands of the in- 
cident electrons: ni = 1 (solid curve), 2 (dashed), 3 (short 
dash), 4 (dash-dot), 5 (dotted). The horizontal axis is the en- 
ergy of the incoming electron measured in terms of the first- 
subband energy of the input lead, E\ = ft 2 7r 2 /2m* ui 2 , from 
a zero of energy at £\. Lower panel: The cooling parame- 
ter n for the coefficients in (a) for initial chemical potential 
fii = 0. The "reduced initial temperature" (horizontal axis) 
is the dimensionless quantity ksTijEx. 



sidearm in Fig. ^ For each geometry we determine the 
initial temperature Ti and chemical potential fii that 
minimize the cooling parameter rj. Table [I] shows these 
data and the corresponding final temperature T a and 
chemical potential /i Q for maximum cooling. 

To illustrate these data and the effect changing the ge- 
ometry in this way, we show in Fig [3] the variation of 77 
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heating (77 > 1) for others. Only the transmission coef- 
ficients depend on the device geometry, so it is through 
these quantum-mechanical scattering probabilities that 
we can control the extent to which a device can cool 
electrons. 



Reduced initial temperature 

FIG. 3: The cooling parameter 77 as a function of the reduced 
initial temperature with initial chemical potential \n — 
for different values of w s /wi. The cooling parameter de- 
creases with increasing magnitude of w s /u>i = 1.0 (solid 
line), 0.9 (long dash), 0.75 (medium dash), 0.6 (short dash), 
0.5 (dash-dot), and 0.4 (dotted). The reduced initial temper- 
ature is defined in the caption to Fig. [2] 
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TABLE I: Cooling parameters n of T-junctions with differ- 
ent sidearm widths w a . Also shown are the input chemical 
potential /ij, output chemical potential fi , and output tem- 
perature To for maximum cooling. 



with initial temperature Ti (for initial chemical potential 
/ij = 0) for different sidearm widths. These results show 
that we achieve cooling (77 < 1) for some geometries and 
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FIG. 4: State-to-lead transmission coefficients for electrons in 
a T-junction device with Wi = w and w s = OAwt. Electron 
are scattered into the output lead from subbands of the input 
lead m = 1 (solid curve), 2 (long dash), 3 (medium dash) 
4 (short dash), and 5 (dotted). The horizontal axis is the 
energy of the incoming electron measured in terms of the first- 
subband energy of the input lead, £\ = ft 2 7r 2 /2m* wf, from 
a zero of energy at £\. See also the data in Tbl. HJ which 
analyzes cooling for this case. 



We shall now consider in detail a device with w s /wi — 
0.4, which gives ~ 0.05% cooling. Figure 0] shows state- 
to-lead transmission coefficients for such a device. Dif- 
ferent curves correspond to different initial subbands. At 
low energies the transmission probability is nearly unity, 
so for this geometry no low-energy electrons are lost from 
the initial distribution. Were we to adjust the Fermi en- 
ergy of this device so only the lowest two subbands were 
occupied, we would see cooling. 

In Fig.[5]we illustrate the dependence of the cooling ef- 
fect on the initial chemical potential p^. This figure shows 
the cooling parameter r\ for w s /wi = 0.4 as a function 
of the initial temperature Ti for /i; = 0.0, 3.0, and 6.0. 
Electron cooling is maximized for /i, = 0, the edge of the 
lowest subband. At this chemical potential all electrons 
scattered into the sidearm are in the thermally active re- 
gion of the Fermi distribution. Cooling is also obtained 
for fii = 3.0, the edge of the second subband. But at 
larger values of \Xi the device heats electrons. 

This example shows that a T-junction can cool elec- 
trons. We would prefer, however, a device that pro- 
duces more cooling than 0.05%. Investigation of other 
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FIG. 5: The cooling parameter ?jasa function of the reduced 
initial temperature for a T-junction device with Wi = w and 
w s = 0.4«7 o for three initial chemical potentials, fii = 0.0 
(solid curve), 3.0 (dashed), and 6.0 (dotted). These data are 
based on the transmission coefficients shown in Fig. 2] Max- 
imum cooling is obtained for fii — 0. The reduced initial 
temperature is denned in the caption to Fig. [5] 

T-junction geometries showed that such a device can- 
not produce significantly more cooling for any value 
of w s /wi. So we next investigate the addition of a sec- 
ond sidearm. This change produces the "plus junction" 
illustrated in Fig. [fj] 

B. A "plus-junction" cooling device 

Since we achieved cooling in a T-junction device with 
w s = 0.4, we shall consider a plus junction with the 
same lead ratios: Wi = w a and w s /w — 0.4. (Note 
that the widths of the two sidearms in Fig. are the 
same.) State-to-lead transmission coefficients for this de- 
vice are shown in Fig. Ufa), and the resulting cooling 
parameter r\ as a function of initial temperature Tj in 
Fig. Efb). The latter figure shows that this geometry 
yields appreciably more electron cooling — greater than 
15% — than did the T-junction devices of Tbl.[I] A single 
unit of a plus-junction device with this geometry would 
cool room-temperature electrons by about AT = — 45C. 
Cooling could be further enhanced by combining several 
such units in sequence. Comparison of this plus-junction 
to the T-junctions discussion previously, as illustrated 
in Fig. [SI show how greatly adding a sidearm improves 
the effectiveness with which this simple device cools elec- 
trons. 

One could try to further optimize the geometry of 



this device by, for example, increasing the width of the 
sidearms with increasing distance from the junction. Al- 
ternatively, one could round the sharp corners at each 
junction into smooth curves. However, our initial explo- 
rations of such alterations (for a T-junction) did not pro- 
duce substantially more cooling than we obtained with 
the far simpler plus junction in Fig. The essential 
features that makes the plus junction more effective than 
any T junction is the presence of more than one sidearm. 
It is also essential the threshold for scattering into these 
sidearms lie above the threshold for the second subband 
of the input lead. 
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FIG. 6: Schematic of a plus-junction. Note that while the 
widths of the input and output leads are equal (wi = w ), 
the width w s of the two sidearms may be smaller or larger 
than Wi. This geometry gives substantially more cooling than 
any device we have explored that has only a single sidearm. 



IV. EXPERIMENTAL CONSIDERATIONS 

We now consider a plus-junction device using realistic 
experimental parameters for Insb and GaAs. While in 
Sec. IIIII we reported results in dimensionless units, we 
here use dimensional units. 

All material properties depend on the Fermi energy Ep 
of the system. This quantity is inversely proportional to 
the effective mass m* of the electrons in the material and 
is determined by the electron density in the reservoir. 
For our device we chose for the initial chemical potential 
(which is approximately equal to the Fermi energy) /ij — 
0. Since we have scaled the energy by S\ and chosen £[ as 
the zero of energy, setting the initial chemical potential 
equals zero means that 

m = E F = 2 , (16) 

2m* wf t 

where Wi is the width of the input lead of the device 
in Fig. The condition l|16fl allows us to determine 
the width of the quantum wire as a function of electron 
density. Since the Fermi energy and the subband energies 
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FIG. 7: Upper panel: state-to-output-lead transmission co- 
efficients for electrons in a plus-junction device with w s — 
OAwi. The curves correspond to different subbands of the 
incident electrons: rii — 1 (solid curve), 2 (dashed), 3 (short 
dash), 4 (dash-dot), 5 (dotted). The horizontal axis is the 
energy of the incoming electron measured in terms of the 
first-subband energy of the input lead, E\ = h 2 ir 2 /2m* w 2 , 
from a zero of energy at £{. Lower panel: Cooling parameter 
r\ based on the coefficients shown in the upper panel. The re- 
duced initial temperature is defined in the caption to Fig. [5] 



depend on the effective mass in the same fashion, the 
width of the quantum wire is independent of the material. 
For a sample with electron density n, we have 



irh 2 



h 2 ir 2 
2m* w 2 ' 



(17) 



FIG. 8: Cooling parameters for a T-junction with u>i = w s = 
Wo (solid line), a T-junction with Wi — w and w„/wo ~ 0.4 
(dashed), and a plus-junction with Wi — w and w 3 = 0.4wi 
(short dash). In all cases the initial chemical potential is \n = 
0. The reduced initial temperature is defined in the caption 
to Fig. H 



which gives for the width of our quantum wire 



Wi 



(18) 



If, for example, 



1.0 x 10 11 cm 2 , we obtain W; 



39.6 nm, quite a small value. We can increase this value 
by decreasing the electron density. 

For a plus-junction, we were able to maximize cool- 
ing by setting the initial temperature to Ti ~ 2 in di- 
mensionless units. In dimensional units, this optimum 
temperature is 



k B T — Topt £{, 
Using Eq. ifTfljl . we obtain 



k B T — T opt Ep. 



(19a) 



(19b) 



For a sample with n = 1.0 x 10 11 cm~ 2 , the initial 
temperature for maximum cooling is Tj ~ 82K for GaAs 
and Ti ~ 399K for InSb. Room temperature (300K) 
corresponds to T ~ 1.5 for InSb. One can therefore ob- 
tain substantial cooling due to quantum effects in a room- 
temperature device. A cooling parameter of r\ ~ 0.9 
implies that the electron population is cooled by 30isT. 
The resulting cooled electrons could be used for photo- 
detection of optical frequencies corresponding to thermal 
energies near room temperature. 

The devices we have considered have only one cooling 
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stage. One could increase cooling by connecting multiple 
plus junctions in series. The spacing between junctions, 
however, must be large enough that resonances in scat- 
tering between sidearms are negligible. If not, one would 
have to treat the device as a single large quantum me- 
chanical scattering target. While the presence of such 
resonances would not preclude cooling, it would make 
calculations for a chain of junction devices more difficult 
and sensitive to details of phase breaking. 



an analytic solution for the amplitudes of outgoing states 
in the leads when only one incoming state is occupied. 




Lead 1 



Lead N 



V. CONCLUSIONS AND PROSPECTS FOR 
FUTURE RESEARCH 



Lead 2 



Many photo-detection applications require a cold de- 
tector. We have presented results for a prototype device 
that demonstrates electron cooling in a single-particle 
picture. We have shown that, while a naive T-junction 
can produce modest cooling and may produce heat- 
ing, adding an additional sidearm yields a device that 
can produce appreciable cooling — at least 15%. The 
abrupt discontinuities in the confining potentials in these 
models are not essential to cooling; what is essential is 
that higher-subband states, which consist of states with 
larger transverse momenta, scatter appreciable into the 
sidearms. We therefore expect electron cooling in such 
devices to be insensitive to details of the potential so 
long as the potential does not eliminate the states of the 
lowest subband. 
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APPENDIX A: 



R-MATRIX THEORY FOR A 2-D 
SYSTEM 



We consider the two-dimensional system in Fig. [5] 
This system has a central region A connected to N exter- 
nal regions or "leads." The leads and the interior region 
meet at a set of boundary surfaces we denote by So, Si, 
. . . , Sn- We treat the boundaries between the shaded 
and unshaded regions as "hard walls" (infinite potential) 
so electron wave functions are non-zero only in the shaded 
regions. Since there may be more than three leads, we 
depart from the notation used in the body of this pa- 
per (in which the input, output and sidearm leads were 
denoted by subscripts i, o, and s) and denote the input 
lead by a zero subscript and all other leads by positive 
integer subscripts. We measure all distances in units of 
wq and energies in terms of Eq = ft 2 /2m* Wq. We seek 



FIG. 9: Schematic of a two-dimensional device for the present 
scattering calculations. The surfaces Si, S2, ... Sn separate 
the interior region A from the N leads. 

The time-independent Schroedinger equation for the 
scattering function is 



H-E)\V 



■E,n 



= 0. 



(Al) 



where \^E,n ) represents the state of an electron with 
kinetic energy E incident in input-lead subband no. Note 
that \^E,n ) is well-defined in all leads. In a, finite region, 
the Hamiltonian H is not Hcrmitian. We can produce a 
Hermitian operator by adding to H the so-called Bloch 
operator Lb~ We denote the eigenfunctions of the sum 
of these operators in the interior region by \<f>i) and write 
the so-called Bloch eigenvalue equation as 



(h + L b ) = Eilfa) 



(A2) 



Inserting the Bloch operator into the Schroedinger equa- 
tion we get 



H + L B -E) * B = L B ^ E . 



(A3) 



We now expand the scattering wave function \^e) hi 
the set of orthonormal Bloch eigenfunctions 



(A4) 



Inserting this expansion into the Schroedinger equation 
and using the properties of the Bloch eigenfunctions 
yields 



I* 



E,n J 



3 



t>j\ L B \*E, n ) , , . 

\<Pj) , 



Ei -E 



(A5) 



where Ej is the eigenvalue corresponds to the Bloch 
eigenfunction \<j)j). This expansion is valid throughout 
the interior region A and on its surface (see Fig. |SJ ■ 
To derive an equation for the R matrix, we now apply 
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this expansion of the scattering state on each boundary 
Si . At each such boundary we can expand the scattering 
function in either lead eigenfunctions or Bloch eigenfunc- 
tions in the interior region. To be specific, we introduce 
a local Cartesian coordinate system for each lead: x q and 
y q are the longitudinal and transverse coordinates of the 
q th lead, respectively. We choose x q = on each bound- 
ary. (One can easily choose any orthonormal coordinate 
system, mutatis mutandis). Each lead eigenfunction is 
then a product of a plane wave in the x q direction and a 
transverse bound-state eigenfunction Xn(Uq)- The scat- 
tering wave function in the q th lead therefore becomes 

^E.n Op, Vp) = e~ ik °' n ° x ° xo,n (yo) S P ,o 

N 

+ T q,n q (Ey kq ^ Xq Xq.n q (y q )6 p , q , 
q,n q = l 

(A6) 

where k q<n and r q >n are the wave vector and transmis- 
sion amplitude for the channel with quantum number n q 
in channel q. Also, Xn q {yq) is the n q h transverse eigen- 
function of lead q. Finally, 8 p _ q is the Kroniker delta- 
function, which ensures that each wave function is de- 
fined only in one lead. If we measure energy in units of 



Eq then we can express energy conservation in lead q as 
E = k 2 q n + n q Tr 2 /w q , where w q is the width of the q th 
lead (in units of wq). We use this equation to determine 
the wave vector k q _ nq . 

After some algebra we get a set of linear algebraic equa- 
tions that we can solve for the transmission amplitudes: 



i^Z T P,n v {E) k p , np Mq >nq ^ np (E) - T q .n q {E) 
P,n p 

+ iko, no M q 

,n q ,0,n ■ 

In writing these equations we have defined 



(A7) 



M q 

,n q ,p,n p — / / Xq,n Q {y q )RE{yq,y P )xp,n p {y P )dy q dy p . 
Jy P Jy q 

(A8) 

Finally, the R-matrix is given by 
R[E,y p ,y q ) = E--E ' ^ ' 

3 3 

This equation is general in that we can easily adapt it to 
any number of leads and to different choices of input lead. 
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